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We demonstrate the quantum fidelity approach for exploring and mapping out quantum phases. As a simple 
model exhibiting a number of distinct quantum phases, we consider the alternating-bond Ising chain using the 
infinite time evolving block decimation method in the infinite matrix product state representation. Examining 
the quantum fidelity with an arbitrary reference state in the whole range of the interaction parameters leads to the 
explicit detection of the doubly degenerate groundstates, indicating a Zt broken symmetry. The discontinuities 
of the fidelity indicate a first-order quantum phase transition between the four ordered phases. In order to char¬ 
acterize each phase, based on the spin configurations from the spin correlations, even and odd antiferromagnetic 
order parameters are introduced. The four defined local order parameters are shown to characterize each phase 
and to exhibit first-order quantum phase transitions between the ordered phases. 
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I. INTRODUCTION 

Tensor network representations have enabled significant 
progress in the computational study of quantum phase tran¬ 
sitions JMl- More specifically, a wave function represented 
by a tensor network is convenient for the simulation of quan¬ 
tum many-body systems. In one-dimensional spin systems, a 
wave function for infinite-size lattices can be described by the 
infinite matrix product state (iMPS) representation 0S]- The 
iMPS and tensor networks in general offer an understanding of 
critical phenomena in infinite and finite lattice systems from 
the perspective of quantum entanglement (see, e.g.. Refs [§- 
[Toil for reviews). It has also been demonstrated that quantum 
fidelity is a useful tool to detect phase transition points and 
degenerate groundstates EHH originating from a sponta¬ 
neous symmetry breaking for a broken symmetry phase, with¬ 
out knowing what type of internal order is present in quan¬ 
tum many-body states. Quantum fidelity based approaches 
have been successfully implemented to investigate quantum 
phase transitions in a number of models. Examples in one- 
dimension include the Ising model in a transverse magnetic 
field 0, the XYX model in an external magnetic field fl3l . 
the bond alternating spin-1/2 Heisenberg chain [15j, and the 
g-state Potts quantum chain |fl3il . 

In this study we further demonstrate the quantum fidelity 
approach for exploring and mapping out quantum phases. As 
a simple but illustrative model exhibiting a number of distinct 
quantum phases, we consider the Ising chain with alternat¬ 
ing interaction strengths. This bond alternating model, also 
known as the dimerized Ising chain, has been used, for ex¬ 
ample, to inve stig ate non-equilibrium spin dynamics at finite 
temperatures lfl6l [l7ll . the so called Glauber dynamics fl8| . 
More recently it has also been investigated in the context of 
additional Dzyaloshinskii-Moriyainteractions Jl9l - i2lll . How¬ 
ever, it has not been fully considered how Landau’s sponta¬ 


neous symmetry breaking picture applies to this model. Here 
we address this issue using the quantum fidelity approach and 
demonstrate how to define the explicit order parameters quan¬ 
tifying the distinct quantum phases throughout the whole pa¬ 
rameter range. 

We calculate the groundstate wavefunction of the infinite 
spin-1/2 bond-alternating Ising chain by employ ing the in¬ 
finite matrix product state (iMPS) representation [|4j, |5|] with 
the infinite time evolving block decimation (iTEBD) method 
developed by Vidal & In order to capture the symmetry- 
broken phases, we use period four matrix product states in¬ 
cluding characteristic eight tensors for the iMPS representa¬ 
tion. The basic idea outlined in Sec. II is to use the quantum 
fidelity with an arbitrary reference state, allowing the doubly 
degenerate groundstates to be detected for the whole range of 
the two exchange interaction parameters. From the disconti¬ 
nuities of the quantum fidelity, we find that a first-order phase 
transition occurs between the ordered phases once one of the 
two interaction strengths changes sign. Further, from the spin 
correlation functions calculated from the degenerate ground 
states, we discuss in Sec. Ill the characteristic spin configu¬ 
ration for each phase and define the possible local order pa¬ 
rameters in Sec. IV, including even and odd antiferromagnetic 
ordering. It is shown that the four defined local order param¬ 
eters reveal four ordered phases and exhibit first-order transi¬ 
tions between the ordered phases in agreement with the results 
using the quantum fidelity. Concluding remarks are given in 
Sec. V. 


II. BOND-ALTERNATING ISING CHAIN AND QUANTUM 
FIDELITY PER SITE 

We consider the spin-1/2 bond-alternating Ising chain 
given by the Hamiltonian 


H = Z{ J ' S 2i-l S 2i + JS 2i S 2 i+ ^ 
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Random Initial State Trials 


FIG. 1: (Color online) Fidelity per site dPV tn) (9), \<p)) between the n- 
th calculated groundstate | V P < " ) (0)) from the n-th random initial state 
and an arbitrary chosen reference state |0) for (a) 9 = n/A, (b) 9 = 
3tt/4, (c) 9 = 5tt/ 4 and (d) 9 = In/A. The horizontal axis denotes 
the number n of the random initial state. It is clearly shown that for 
each case there exist two degenerate groundstates corresponding to 
the two values of d. 


where St is the spin operator on the z'th site. The exchange 
couplings are J and J', which we parametrise in terms of the 
variable 9, with J = cos 6 and J' = sin 9. For J — J' < 0 
(9 = 5n/4) the system becomes the ferromagnetic (FM) Ising 
model and for J - J' > 0 (9 - n/4) the antiferromagnetic 
(AFM) Ising model. In the both cases J - J' > 0 and 
J — J' < 0 the Hamiltonian is one-site translational invariant. 
Due to the bond alternation, then, the Hamiltonian is two-site 
translational invariant, except for the cases J - J' > 0 and 
J - J' < 0. 

An iMPS ground state of the system can be obtained by 
using the iTEBD algorithm with a chosen initial state in 
the iMPS representation. We calculate n groundstate values 
| '¥ (n \9)) corresponding to the / 7 -th random initial state. In or¬ 
der to determine how many groundstates exists for a given pa¬ 
rameter 9, we consider the quantum fidelity F(\ x ¥ (n \9)), \<p )) = 
|('p(”*(0)|0)| which is the overlap function between the n- 
th calculated ground state | V F (,,) ( 0 )) and an arbitrary refer¬ 
ence state | (j)). For our numerical study, the reference state 
| <p) is chosen randomly. The quantum fidelity scales as 
F (|ip«( 0 )\i 0 » ~ d L , where L is the system size. Following, 
e.g.. Ref. [Tj, the fidelity per site can be defined as 

im (2) 

L—> oo 

From the fidelity F('P ( " ) (0)), |</>)), the fidelity per site satisfies 
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FIG. 2: (Color online) Fidelity per site between the doubly de¬ 
generate groundstates with the same reference state used in Fig. |T| 
as a function of the interaction parameter 9. The four different or¬ 
dered phases {I. II, III, IV} are clearly distinguished. For each ordered 
phase, the two values of the fidelity per site indicate a doubly degen¬ 
erate groundstate. 


the properties: (i) normalization d{\<p), |</>)) = 1 and (ii) range 

0 < r/(T' ( " ) (60>, |0» < 1. 

A degenerate groundstate can be determined from the fi¬ 
delity per site r/(| v P ( ' ,) ), |</>)) as a function of the random initial 
state trials n, as shown in Fig.Q]for (a) 9 = n/4, (b) 9 = 3n/4, 
(c) 9 = 5n/4, and (d) 9 - ln/4. For the iMPS representa¬ 
tion, the truncation dimension^ used is x - 32. These plots 
show that there are two different values of the fidelity per site 
for the groundstates from the 30 random initial states. For a 
large enough number of random initial state trials, the prob¬ 
ability P(n) of each degenerate groundstate approaches 1/2, 
i.e., lim^oo P(n) = 1/2. This implies that there are doubly 
degenerate groundstates for each given 9 value. 

In order to determine how many ordered phases there are in 
the model, we calculated the fidelity per site as a function of 
the interaction parameter 0 < 9 < 2n (see Fig. [2}. Here, we 
have chosen the same reference state |</>) as in Fig. Q] Fig. [2] 
shows clearly that the fidelity is discontinuous at four points, 
i.e., 0 = 0, 7 t/ 2, 7r, and 3nj2. The discontinuous fidelities 
indicate that a first-order quantum phase transition occurs at 
each discontinuous point corresponding to each critical point. 
Consequently, the system has four ordered phases: 

I: 0 <9 < nil (J > 0 and J' > 0), 

II: n/2 < 9 < n (J < 0 and J' > 0), 

III: n < 9 < 3n/2 {J < 0 and J' < 0), 

IV: 3n/2 <9<2n(J>Q and /' < 0). 

Each phase has doubly degenerate groundstates denoted by 
| x H Q .,j(6»)> and |T, K , 2 ((0)), where a e {I, II, III, IV} labels the 
phases. According to the spontaneous symmetry breaking the¬ 
ory, the doubly degenerate groundstates imply that a Z 2 sym¬ 
metry is broken for each phase. Hence, for each phase, a dif¬ 
ferent Z 2 symmetry is broken and the system state belongs to 
a different ordered phase. 
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FIG. 3: The spin correlation O z ss (i, j) as a function of the site distance 
|i - j\ for (a) 9 = n/ 4, (b) 9 = 3n/A, (c) 9 = 5tt/ 4, (d) 9 = In /4. 


III. SPIN CORRELATIONS AND GROUNDSTATE 
WAVEFUNCTIONS 


In order to gain further insight into the four ordered phases, 
we consider the spin-spin correlations for each phase. The 
spin-spin correlation function is defined by 


written as 


oo 


1^1,1(21/4)) = Y] IT>2i Il>2, + 1, 

i ——00 

(4a) 

00 

1^1,2(21/4)) = Y\ lt»2f+l • 

(4b) 


1=1 — 00 


Thus, for 0 < 6 < n/2 (J > 0 and /' > 0), the groundstate 
is in the antiferromagnetic (AFM) phase. The groundstate 
wavefunctions are two-site translational invariant. However, 
the Hamiltonian in Eq. (|T|) is two-site translational invariant 
for / 4 J' but is one-site translation invariant for J - J'. 
Also, note that one of the two groundstates transforms to the 
other groundstate under the spin-flip transformation. Hence, 
for J 4 J', the two groundstates result from the spontaneous 
symmetry breaking of the spin-flip symmetry. We can think of 
the two groundstates in the AFM phase as losing more sym¬ 
metry for J - J’ than J 4 J' from the spontaneous symmetry 
breaking of the spin-flip symmetry. 


B. Odd antiferromagnetic phase 

For 6 = 3n/4, the spin correlation in Fig. [3jb) implies that 
the spin configuration is either T/T;+il;+ 2 l;+ 3 T ;+4 ••■or 
• • • I 1 I/+ 1 T 1 + 2 T/+ 3 I /+4 • • • ■ The groundstate wavefunctions are 

OO 

= Y] lt>4i IT>4i + l ll>4,+2 ll>4, + 3 . (5a) 

i ——00 
00 

I^H2(3w/4)> = Y] U>4i Ii>4« + 1 lt>4/+2 lt>4/-3 • (5b) 


O z ss (i, j) = (07 cH), (3) 

where i and j denote locations along the chain. In Fig. [3j 
we plot the spin-spin correlation O z s (i, j) as a function of the 
separation distance |i - j | for the typical parameter values (a) 
6 — 7t/4, (b) 9 = 3n/4, (c) 9 - 5n/4, and (d) 9 — ln/4 in each 
of the phases. 


A. Antiferromagnetic phase 

The two groundstates for each phase give the same spin- 
spin correlation. Furthermore, the spin-spin correlation as 
a function of the lattice distance has only two values, i.e., 
+ 1 or -1. This result implies that the direction of the two 
spins i and j is either parallel for 0 7 ss (\i - j\ ) = +1 or anti¬ 
parallel for O z s (\i - /I) = -1. For instance, for 9 = 7r/4, 
the spin configurations, using an obvious notation, are either 
T/X/+ 1 T 1 + 2 Xi+ 3 Ti +4 or ■ • • XzT i+ 1 X/+ 2 Tz+ 3 Xz +4 •• Con¬ 

sequently, for 9 = n/4, the groundstate wavefunctions can be 


Thus, for n/2 < 9 < n (J < 0 and /' > 0), the ground- 
state wavefunctions are four-site translational invariant be¬ 
cause the exchange interaction on the even bonds is ferro¬ 
magnetic (FM), i.e., / < 0, while the exchange interaction 
on the odd bonds is AFM, i.e., J' > 0. Note that one of 
the two groundstates transforms to the other groundstate un¬ 
der the spin-flip transformation and the two-site translational 
transformation. The Hamiltonian is invariant for the spin-flip 
transformation and the two-site translational transformation. 
Hence, the two degenerate ground states arise from the spon¬ 
taneous symmetry breaking of the group G — ZiX Z 2 /Z 2 , i.e., 
the spin-flip Z 2 and the two-site translational Z 2 symmetries. 
The spin configuration shows two-spin alternating behavior. 
In order to distinguish this phase from the AFM phase in re¬ 
gion I, we then call this phase the odd AFM phase. 


C. Ferromagnetic phase 

The spin correlation for 9 = 5n/4 in Fig. [^c) isFM. For 7 r < 
9 < 3n/2 (J < 0 and J' < 0) the groundstate wavefunctions 
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are the FM wavefunctions 


oo 


l'Fm,i(57r/4)> = Y] IT>, 

i=— 00 

(6a) 

00 

l'Pm,2(57r/4)) = f] U), . 

(6b) 


l——CO 


The groundstate wavefunctions are one-site translational in¬ 
variant. However, the Hamiltonian in Eq. ([T]) is two-site trans¬ 
lational invariant for / + J' but is one-site translation invariant 
for ] - . Also, note that one of the two groundstates trans¬ 

forms to the other groundstate under the spin-flip transforma¬ 
tion. Hence, for J 4 J' , the two groundstates result from 
the spontaneous symmetry breaking of the spin-flip symme¬ 
try. For J = interestingly, the two groundstates in the FM 
phase have more symmetry than the Hamiltonian because the 
two groundstates have one-site translational symmetry but the 
Hamiltonian has two-site translational symmetry. Such a sym¬ 
metry, which is not directly manifest in the Hamiltonian, is an 
example of an enhanced or emergent symmetry |[22l - l28i] . 

D. Even antiferromagnetic phase 

The spin correlation for 6 = In /4 in Fig. 0d) is similar to 
the case 6 = 3n/4. Here the groundstate wavefunctions are 

OO 

mv,i(7TT/4)> = n ^11)4,1,11)4,2 lt)4, + 3, (7a) 

i=—oo 

OO 

l^iv, 2 (77r/4)) = ( 7b > 

/=—OO 

Thus for 37 t /2 < 9 < 2n (J > 0 and J' < 0) the ex¬ 
change interaction on the odd bonds is FM ( J ' < 0) while 
the exchange interaction on the even bonds is AFM (J > 0). 
The groundstate wavefunctions are four-site translational in¬ 
variant. Note that one of the two groundstates transforms to 
the other groundstate under the spin-flip transformation and 
the two-site translational transformation. Also, the Hamilto¬ 
nian is invariant for the spin-flip transformation and the two- 
site translational transformation. Hence, the two degenerate 
ground states arise from the spontaneous symmetry breaking 
of the group G - Z 2 x Z 2 /Z 2 , i.e., the spin-flip Z 2 and the 
two-site translational Z 2 symmetries. The spin configuration 
shows two-spin alternating behavior. In order to distinguish 
this phase from the AFM phase in the regions I and II, we 
then call this phase the even AFM phase. However, both of 
the odd and even AFM phases originate from the breaking of 
the same symmetry. Then one needs to consider how to distin¬ 
guish the two phases. This point is clarified in the discussion 
of the appropriate order parameters in Sec. IV. 

IV. ORDER PARAMETERS 

To distinguish the four phases, one has to define character¬ 
istic local order parameters. Based on the spin configurations 



0/tc 


FIG. 4: (Color online) The four order parameters as a function of 
the interaction parameter 0. In (a)-(d), the blue triangles and the red 
circles denote the average value of the local order parameters, defined 
in the text, with each ground state, | v P t , > i) and flf,^), respectively. 


or the wavefunctions, we define the order parameters for the 
four phases. Normally, one uses the FM and the AFM order 
parameters defined by the magnetization M F — {a~ + cH +1 )/2 
and the staggered magnetization Mf F — {cr\ - o~ +1 )/2. How¬ 
ever, for the bond-alternating model, these definitions for the 
order parameters do not distinguish between all four phases, 
because of the even AFM and odd AFM phases. Further¬ 
more, the odd and even AFM phase originate from the spon¬ 
taneous breaking of the same symmetries. To overcome this, 
one needs to consider the symmetries of the groundstates. 
One can notice that from the two degenerate groundstates for 
each phase, the defined FM order parameter should be one- 
site translational invariant, the AFM order parameter two-site 
translational invariant, and both the odd and even AFM or¬ 
der parameters four-site translational invariant. Then, proper 
order parameters should be four-site translational invariant. 
From the properties of the groundstates for each phase, in 
terms of the normal definitions of the magnetization and the 
staggered magnetization, we define four order parameters as 
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FIG. 5: (Color online) Phase diagram for the spin-1/2 bond- 
alternating Ising chain. 


follows, 


M AFM ~ 

1 

2 

+ m£ 2 ). 

( 8 a) 

M'fm 

1 

2 

K+< +2 ). 

( 8 b) 

\; ' 

LV1 AFM, even 

1 

2 

«-a4 +2 ), 

( 8 c) 

M Z — 

1V1 AF M,odd 

1 

2 

K +1 

( 8 d) 


We plot these order parameters as a function of the interac¬ 
tion parameter 8 in Fig. [4] It is shown clearly that each defined 
order parameter is non-zero for the region of each of the or¬ 
dered phases, otherwise zero. This shows that the defined or¬ 
der parameters characterize each of the ordered phases. Fur¬ 
thermore, the quantum phase transitions are seen to be first 
order by the sudden drop of the order parameters to zero at 
the discontinuous points of the fidelity revealed in Fig. [2] At 
these points the groundstates are infinitely degenerate. As a 
result, the spin-1/2 bond-alternating Ising model has the four 
ordered phases. As is shown in the phase diagram of Fig. 0 
the system is clearly seen to be in the 

(i) AFM phase for 0 < 8 < n/2 {J > 0 and J' > 0), 

(ii) odd AFM phase for 7r/2 < 8 <n {J < 0 and J' > 0), 

(iii) FM phase for n < 9 < ln/2 (J < 0 and J' < 0), 

(iv) even AFM phase for 37t/2 < 6 < 2n (J > 0 and J' < 0). 

The phase boundaries are the lines J - 0 and J' = 0. 

Actually, if the system undergoes a full spontaneous sym¬ 
metry breaking of the group G = Zi X Z 2 , there are four de¬ 
generate groundstates. However, we have detected only the 


two degenerate groundstates in each phase from the fidelity 
shown in Fig. [2] and the order parameters in Fig. [4] As one 
can notice in Fig. [4] the odd and even AFM phases are clearly 
distinguished by the defined order parameters. These show 
that for the odd and even AFM phase, the system undergoes 
a partial spontaneous symmetry breaking of the group, which 
induces the two degenerate ground states for each phase. De¬ 
pending on the sign of the interaction strengths J and J', the 
partial spontaneous symmetry breaking generates the different 
ordered phases, i.e., the odd and even AFM phases. Conse¬ 
quently, the odd and even AFM phases are distinguishable. 


V. CONCLUSION 


We have investigated quantum fidelity in an infinite-size 
bond-alternating Ising chain by employing the iMPS represen¬ 
tation with the iTEBD method. By detecting the doubly de¬ 
generate groundstates for each phase by means of the quantum 
fidelity with an arbitrary reference state, it was shown that, for 
each phase, a different Z 2 symmetry is broken and the system 
state belongs to a different ordered phase. By also detecting 
the discontinuities of the quantum fidelity, we demonstrated 
that first-order quantum phase transitions occur between the 
ordered phases as the interaction parameter 8 varies through 
0 < 8 < 2n. Based on the spin configurations from the char¬ 
acteristic properties of the spin correlations, the four defined 
local order parameters, including the even and the odd AFM 
order parameters, are shown to clearly characterize each phase 
and the existence of the first-order quantum phase transitions 
between the ordered phases. Consequently, by taking a sim¬ 
ple and well known model as example, we have demonstrated 
the usefulness of the quantum fidelity with an arbitrary refer¬ 
ence state to investigate the nature of quantum phases without 
knowing a priori what type of internal order is present in a 
quantum many-body state. 
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